Swing switching of spin-torque valves 
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We propose a method for inducing magnetization reversal using an AC spin current polarized 
perpendicular to the equilibrium magnetization of the free magnetic layer. We show that the critical 
AC spin current is significantly smaller than the corresponding DC one. The effect is understood 
as a consequence of the underdamped nature of the spin-torque oscillators. It allows to use the 
kinetic inertia to overcome the residual energy barrier, rather than suppressing the latter by a large 
spin current. The effect is similar to a swing which may be set into high amplitude motion by a 
weak near-resonant push. The optimal AC frequency is identified as the upper bifurcation frequency 
of the corresponding driven nonlinear oscillator. Together with fast switching times it makes the 
perpendicular AC method to be the most efficient way to realize spin-torque memory valve. 

PACS numbers: 75.70.-i, 85.75.-d, 75.75.Jn 



I. INTRODUCTION 

Spin-torque (ST) induced switching in magnetic tun- 
nel junctions (MTJ) has been a subject of intense exper- 
imental and theoretical research since its inception by 
Slonczewski [l[ and Berger [2j over a decade ago. It has 
been demonstrated [3 6] to be a viable candidate for fast, 
scalable and non-volatile memory storage. Early MTJ 
devices 04111 utilized two ferromagnetic layers, sepa- 
rated by a non-magnetic layer, with magnetization axes 
aligned parallel to each other and orthogonal to the verti- 
cal axis of the pillar, Fig. []Ji (hereafter referred to as the 
parallel method). These devices tend to have switching 
times longer than Ins [l^ - [l4l ] as well as a broad switching 
time distribution [15j . These characteristics come from 
an initial incubation time during the switching process 
when thermal noise provides for an initial incipient mis- 
alignment between the two layers 1161 Il7| . Sub nanosec- 
ond switching times are possible [18j . but require rather 
strong currents and are far from the optimal efficiency 
[l9j . It was recently proposed by Kent et all [2(| that 
adding a third ferromagnetic layer, with an easy axis per- 
pendicular to the other two, i.e. along the vertical axis of 
the MTJ, Fig. [TJ), may reduce switching times below Ins 
[13l [lj, [2l|, [22| ■ The improvement is due to the fact that 
the torque acts perpendicular to the initial magnetization 
direction and therefore does not rely on thermal noise to 
initiate the switch. This approach, here known as the DC 
perpendicular method, allows magnetic reversals of the 
free layer at a fraction of the energy cost, i.e. associated 
Joule heating 22]. While a significant improvement over 
the parallel method, it has it's own deficiency. The rea- 
son is that the magnetization tends to precess around the 
easy axis direction. As a result the perpendicular torque 
acts in the wrong direction for a half of the precession pe- 
riod [2(| • Therefore the protocol is most efficient if it can 
accomplish the switch during the half of the precession 
period. This dictates rather high values of the required 
spin current. 

In this paper we suggest an alternative approach to 
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FIG. 1: (Color online) (a) Diagram of parallel MTJ. (b) Di- 
agram of perpendicular MTJ. Blue and red regions represent 
fixed and free layers respectively. Yellow regions are non- 
magnetic spacers. Current flows along the vertical axis of the 
pillars. 



magnetic reversal in MTJs using an AC current in the 
perpendicular configuration (here known as the AC per- 
pendicular method) with the driving frequency near the 
natural precession frequency of the free layer. The use 
of microwave frequency AC spin-current to influence the 
magnetization dynamics of the free layer has been studied 
in a number of works in recent years. References 0, [24[ 
employed a weak AC signal to excite ferromagnetic res- 
onance excitations of the free layer, detected through 
the induced DC voltage in the perpendicular geometry. 
References (25l . [26j j performed similar experiments in the 
parallel geometry. References [13, HH showed that in the 
parallel geometry an AC signal on top of a DC pulse may 
facilitate stochastic thermal switching events between the 
two metastable states. And reference .29] showed that an 
AC signal prior to a DC pulse can also facilitate switching 
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by reducing the critical current in the parallel geometry. 

Here we show that a properly tuned purely AC current 
in the perpendicular geometry may serve as most energy 
efficient way to switch ST memory valves. Our simu- 
lations show that this method produces fast switching 
times, below Ins, and allows for switching currents about 
a factor of three below the critical current of the equiv- 
alent DC device. As a result it accomplishes the switch 
with a significantly smaller energy dissipation than the 
DC method. We develop a theory which explains the be- 
havior observed in the simulations and serves as a useful 
guide to optimize the system parameters. It is based on 
the fact that the energy of the Stoner-Wohlfarth (SW) 
orbital and the relative phase between the magnetization 
and the external AC drive form a canonical action-angle 
pair [30]. As a result, the system may be described as a 
weakly damped non-linear oscillator. Its effective poten- 
tial landscape is affected by both the amplitude and the 
frequency of the AC spin current. The key observation, 
explaining the critical current reduction, is the inertial, 
i.e. underdamped, nature of this effective oscillator. It 
allows the switching to occur without fully suppressing 
the energy barrier, but rather overcoming it by inertia. 
The effect is similar to setting a swing into a rotation by 
weak near-resonance AC push. 
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where Tll is the torque induced by an effective magnetic 
field which incorporates both the conservative and ther- 
mal stochastic components. Gilbert damping Tqd phe- 
nomenologically describes the rate with which the energy 
is removed from the system. Its torque is perpendicular 
to the Landau-Lifshitz one as well as to the SW orbits of 
the constant energy. The last term rsi on describes the 
effect of the spin current on the magnetization direction 
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Here I s (t) = e p T s (t) is the time dependent spin current 
in units of magnetization, polarized along the direction 
e p . For the DC perpendicular method I s (t) is a rectan- 
gular pulse with the duration t v \ and the amplitude I s , 
while for AC case we take it as I s sin(wi) . An example of 
the AC spin current pulse is shown in Fig. ©a. For most 
simulations discussed here we took the pulse duration as 
three full periods t p \ = 3(2it/uj) of AC signal. 



II. SIMULATIONS OF THE PERPENDICULAR 
SPIN-TORQUE VALVE 
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FIG. 2: (Color online) Switching probability of DC perpendic- 
ular valve as a function of spin current amplitude (X s /M 3 ) and 
an inverse pulse duration t p \ (in units of -yMs) with H x = Ms, 
H z = 0.033M 3 , a = 0.015, and T = 300A. 



To model the magnetic switching we treat the free layer 
as a single magnetic domain with a constant saturation 
magnetization M s and magnetization direction specified 
by a time-dependent unit vector m(i). Its dynamics is 
described by the Landau-Lifshitz-Gilbert equation with 
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FIG. 3: (Color online) Switching probability of AC perpen- 
dicular valve as a function of spin current amplitude (I s /M s ) 
and AC frequency uj/2h. The pulse duration is taken to be 
three driving periods t p i = 3(2n/u}). Same parameters as in 
Fig.H 



The effective magnetic field is given by 
1 



H eff (m,i) 
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Here E(rh) is the conservative anisotropy energy, which 
we take as 
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where H z and H x are the easy and hard axis anisotropy 
fields respectively. The energy is zero E = along the 
easy axis directions ±e z and reaches maximum E = 
fioM s (H z + H x )/2 along the hard axis directions ±e x . 
There are two saddle points along the ±e y axis, which 
separate the two minima and act as the energy barrier 
with the energy Eb = fioM s H z /2. The magnetization 
must overcome this barrier before the switching occurs. 
Thermal noise is included via a random Gaussian mag- 
netic field h(t) [3ljj with the correlator determined by the 
fluctuation-dissipation theorem [34| 



with longer pulse times. 
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where T is the temperature of the free layer. 

Using the above model we performed numeric simu- 
lations of switching processes for DC-perpendicular and 
AC-perpendicular spin currents. We first present results 
for material parameters chosen as M s = 1200Emu / cm 2 , 
H x = M s , H z = 0.033M S , a = 0.015 and T = 300iC 
Other choices of parameters are discussed after the theo- 
retical part. In both DC and AC perpendicular cases the 
spin current is polarized along the z and x axes e p = x+z 
to simulate the influence of the two fixed layers, Fig. [1}d. 

For the DC case Fig.[2]shows the switching probability 
as a function of spin current amplitude T s and the inverse 
pulse duration l/t p \. One may notice that the switching 
current saturates at about I s = 0.014M S in the limit of 
long pulses (the inverse period of small energy precessions 
is about 0.037-Ms). As was mentioned in the introduc- 
tion the switch is most reliable then the pulse duration 
is less than half of the precession period. In this case the 
perpendicular torque is always acting in the proper direc- 
tion of increasing the precession amplitude. For shorter 
pulses the critical current increases roughly as t~f, the 
fidelity of the switch improves as well. 

For the AC-perpendicular case Fig.[3]shows the switch- 
ing probability as a function of spin current amplitude I s 
and AC frequency uj/2it with the pulse duration of three 
oscillation periods t p \ — 3(2tt/lu), see also Fig. [5^. No- 
tice that at optimal conditions the switching occurs at 
the current nearly three times lower than that in the DC 
perpendicular case with the switching currents as low as 
I s = 0.005Af s . One should also notice the sharp fre- 
quency selectivity of the switching process. The optimal 
AC frequency appears to be somewhat below the natural 
small precession frequency £Iq = (27r) x 0.037M S . 

Simulations were also performed for pulses with dif- 
ferent number of spin-current oscillations t p \ — n(2-K/ui). 
For pulses with n < 3 regions of low current switching 
were gradually excluded as n approached 1/2 until the 
switching probability was qualitatively the same as in 
the DC case of Fig. [3] For pulses with n > 3 the range 
of currents and frequencies which produced switching re- 
mained qualitatively the same at in Fig. [3J However, 
the switching probabilities are reduced and oscillate with 
increasing n reaching probability of 50% for sufficiently 
large values of n. This behavior is the result of magneti- 
zation switching back and forth between the two minima 
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FIG. 4: (Color online) Simulated optimal Joule heat lost 
during switching process in arbitrary units as a function of 
spin current amplitude (X s /M a ) for DC perpendicular (black 
squares) and AC perpendicular (red triangles) spin current 
methods. For AC simulations frequencies u)/2tc ranged from 
— 0.037-Ms with pulse times t p i = n(2n/uj) where n ranged 
from 0.5 — 4 to find the optimal regime. Same parameters as 
in Fig. H 

Fig. [1] shows the minimal Joule heat lost during the 
switch as a function of the spin current magnitude for 
both the DC and AC perpendicular cases. It is cal- 
culated by numerically integrating the power over time 
cx R Jq F ' I 2 (t)dt, where R is a resistance which we as- 
sume to be constant through the switching process. Only 
combinations of spin current T s and u which produced 
switching probabilities above 90% were included. Notice 
the AC perpendicular method cost significantly less en- 
ergy to produce switching at low spin current strengths. 
At it's most energy efficient point the AC perpendicular 
method costs nearly one fourth the energy of the the most 
energy efficient DC perpendicular method. This point 
occurs at driving frequency w/27r ~ 0.0227M S , spin cur- 
rent amplitude T s ~ 0.005M S , and the pulse duration of 
< pl = 3(2tt/w)- 



III. ANALYTICAL CONSIDERATIONS 

We look now for a qualitative understanding of why 
the AC protocol requires a significantly smaller critical 
spin current. Another feature calling for an explana- 
tion is a sharp frequency dependence of the switching 
efficiency seen in Fig. [3J Figure [5^, shows the AC spin 
torque pulse T a (t) and the time dependence of the energy 
E(t) = E(ih(t)) along a switching trajectory. Notice that 
the energy increases up to a bit over the top of the energy 
barrier E\, during the time of the AC pulse and later de- 
creases down to zero in the state with the reversed mag- 
netization. Notice also that the energy E(t) is a rather 
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smooth function of time. On the other hand the mag- 
netization rh(t), depicted in Fig. [5b, goes though three 
full revolutions before the switching, exhibiting highly 
non-monotonic behavior. 
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FIG. 5: (Color online) a) Spin current T s (t) — X„cos(u}t) in 
arbitrary units (red, dashed) and normalized energy E/Ei, 
(black, solid) vs. time (in units of (7-M s ) _ ) for a success- 
ful switching event, b) Magnetization M z /M s (black, solid) 
and M x /M a (red, dashed) in units of (jMs)' 1 vs. time for 
the same event. Here l a = 0.005A'/ S and co/2tv = 0.0227M S . 
Other parameters are the same as in Fig. [2] 

This observation suggests to parameterize the magne- 
tization vector rh in terms of a slow variable E, which is 
nothing more than the energy of the free layer given by 
equation (j4]), and a fast variable ip which runs along a 
closed SW orbit of a constant energy. It is convenient to 
make the latter linearly related to the physical time (in 
the absence of the spin-torque). To this end we identify ip 
with the normalized time required to reach a given point 
along SW orbit within the purely conservative evolution, 
i.e. if drh/dt = Tll (we omit the thermal stochastic 
noise ((S|) throughout this section). With this choice one 
takes dip = fl(E)dt, where, cf. Eq. ©, 
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and 0(.E) is the frequency of SW orbit, given by 
2tt/Q(E) — § dt, where the integral runs along the SW 
orbit with energy E and dt is given by Eq. ([5]). This 
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Notice that the parametrization is locally orthogonal, 
i.e. c^m • c>Erh = 0. One can thus project the iden- 
tity drh/dt — dtpfh ip + (?£;m E onto these two orthogonal 
directions and employ LLG equation ([1]) for drh/dt to 
find equations of motion for E and ip pair 
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Here the right hand sides may be also expressed as func- 
tions of E and ip. The energy evolution is driven only by 
relatively weak damping and spin-torque terms, making 
the energy a slow variable. On the other hand, the phase 
dynamics is primarily given by a large il(E), coming from 
LL term and making the phase the fast variable. 

In presence of near resonant AC spin-torque I s sin(wt) 
the fast rotation phase ip tends to follow the external 
drive, see Fig. [5] One can thus introduce the relative 
phase <p(t) = ip(t) — tot as another slow variable. The 
equation of motion for the slow phase </>(i) follows triv- 
ially from Eq. (0) and the relation 



ip — to 
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One can now employ the fact that both E and (f> change 
relatively little through one precession cycle to average 
over the fast precessing variable ip 0,H^|. To this end we 
write the AC spin current as X s sin(wi) = T s sin(</?— <p) = 
T s (sin ip cos 4> — cos ip sin <fi) and integrate the right hand 
sides of Eqs. ©, ([TO]) over dp with the help of Eq. 
This way we arrive at the pair of equations of motion for 
the slow degrees of freedom 



-aU(E) - X S V{E) sin^, 
Q(E) -uo-X s W{E) cost 
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Here the generalized forces averaged over the orbit with 
energy E are given by 



U(E) 
V(E) 
W{E) 



«/MMxH cff ].M, 

j[dM x M] • e p cos ip . 

n(E) 2 i (dM • r LL ) (r LL • e p ) 



(12) 



7 



2ttM s 



■ LL 



■ sin ip, 



where in the last two expressions we used the parity with 
respect to ip to single out odd and even components. 

A typical parametric trajectory E{t) versus 4>{t) is 
shown in figure ^ for a subcritical spin current. No- 
tice that after averaging over the fast fluctuations E(t) 
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FIG. 6: (Color online) Simulated energy E vs phase cj> (black, 
solid) calculated from numerical integration of the LLG equa- 
tion ([1} and average energy and phase (red, dashed) calcu- 
lated from integration of equations Here X s = 0.004M S 
and uj/2-k = O.O287M.,. Other parameters are the same as in 
Fig.H 



FIG. 7: (Color online) Precession frequency S}(E) = uj (red, 
dashed) along with the equilibrium energy E eq (black, solid) 
vs. driving frequency uj calculated from a) numeric integra- 
tion of the LLG equation, Eq. |T} and b) Eqs. (fl3l) . At 
frequencies uj\ and uii the simulation jump between high and 
low energy driven precession modes. Here T s — 0.002M S and 
T = 0. Other parameters same as in Fig. [21 



and 4>{t) behave similarly to coordinate/momentum pair 
of an underdampcd oscillator. The oscillation frequency 
of this effective oscillator is much smaller than the pre- 
cession frequency Q(E) and is rather given by tt(E) — uj. 
Of particular importance is the fact that before wind- 
ing down into the stationary point (E eq , </>cq), the energy 
highly overshoots its equilibrium value. Had the energy 
exceeded the height of the energy barrier Eb in the course 
of such an overshoot - a switch is likely to occur. This 
remains to be the case even if the equilibrium value E eq 
is still below the barrier hight Ej,. This is the qualita- 
tive reason for the AC critical current being significantly 
smaller than the DC one. 

To make progress in explaining sharp frequency depen- 
dence of the effect, one needs to look closely at the equi- 
librium point (E eq , 4> eq ) of the effective oscillator. Such 
a point would be reached by the free layer magnetization 
upon continuous AC drive with the amplitude below the 
switching one after all transients decay due to dissipa- 
tion. Setting the left hand sides of equations (fTTj) to zero 
gives such equilibrium conditions 
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Figure [7J shows the equilibrium energy E cq over a range 
of frequencies along with simulated result found by adia- 
batically varying w over time. Notice that for the range 
of frequencies uj\ < uj < u>2 there are two stable solu- 
tions and an unstable branch in between. Consequently 
the simulation reveal hysteretic loop with two jump fre- 
quencies, where the equilibrium energy abruptly jumps 



to a higher or lower value. Since the goal is to reach the 
switching energy E^, one wants to insure that the oscil- 
lator follows the upper of the two branches, while having 
as large equilibrium energy as possible. As seen from 
Fig. [7] this is achieved at frequencies right above the up- 
per hysteresis jump u>2 < uj. We thus calculate the upper 
jump frequency of the hysteresis as a function of the spin 
current T s and plot in Fig. [§] 

To find a critical current of the swing switching it is 
useful to notice that in the absence of the dissipation 
a — > 0, the AC spin-torque driven oscillator (fTTj) pos- 
sesses an integral of motion. Indeed, one may check that 
the following function 



dE'J(E') uj - fl(E') + l s W(E') cos < 



(14) 

is conserved by the equations of motion (jlip . with a = 0, 
if the function J{E) is a solution of the following linear 
homogeneous differential equation: 

v«©-flDfw(D-™). (15) 



dE 



In fact one may define the action 1(E) — J dE'J'(E'), 
such that the change of variables (E, </>) — > (I, <j>) in 
Eq. (fTT)) results in the Hamiltonian H(I,(f>) written in 
terms of the canonical action-angle pair. The equations 
of motions fTTj are nothing but the Hamilton equations: 
/ = d^H. and <j> = —di%. 

The consequence of this observation is that without 
the dissipation the E, (ft trajectories are closed. The cor- 
responding phase portraits are are plotted in Fig.[5^i,b for 
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FIG. 8: (Color online) Contours of a constant T-L for X a = 
0.003M S and for driving frequencies a) to < 0J2 (uj/2-k — 
0.027M S ) and b) lj > lj 2 (uj/2tv = 0.028jM s ). Separatrix 
lines H — are shown in bold red. Points A, B, and C corre- 
spond to equilibrium values E eq and (j> eq , cf. Fig. [7] Dashed 
black line in b) is a simulated trajectory for same parameters, 
already shown in Fig. [6] 



the case of uj < u>2 arid u> > UJ2 , correspondingly. There 
are two stable points A and B in the former case and only 
one C in the latter. Hereafter we focus on the latter case, 
since for the former the evolution starting from E = 
can only bring the system to a small energy fixed point 
B and thus is not beneficial for the switching. Starting 
from a vicinity of E = line the system evolves along 
a trajectory which is close to the separatrix line H = 0, 
before winding down towards the fixed point C due to 
dissipation. For the switching to take place one must 
require that the separatrix goes beyond the top of the 
barrier E^. Since the maximum extent of the separatrix 
is reached at <j> = 7r, a necessary condition for switching 
to occur is: 



I, > 



J^dEJ(E)(u-n(E)) 
J Eb dEJ(E)W(E) 



(16) 



In fact this is an underestimate, since it does not take into 
account the dissipation. The latter dictates that the sep- 
aratrix must exceed Ej, + Ed, where Ed = a J dtU(E) = 

a J Eh dEU (E) / E is the work done by the dissipative 
force along the separatrix trajectory. One can now sub- 
stitute E = —T S V{E) sin cf>, where sin is found from the 
separatrix equation H = 0, to find the dissipative loss 
Ed- One then uses Ef, + Ed instead of Ef, in Eq. (fT6|) 
to obtain the minimal switching current. The resulting 
critical line I s {u}) is plotted in Fig. [9] One can see that 
this line is indeed close to the boundary of the observed 
switching range. 

To make the general theory more transparent it is use- 
ful to evaluate the averaged generalized forces (|12[) at 
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FIG. 9: (Color online) Critical frequency cj = u>2 (dashed) and 
critical current (dot-dashed) lines atop of the switching prob- 
ability as a function of spin current (I 3 /M a ) and driving fre- 
quency (in units of "(Ms) for the pulse duration t p \ = 3(27t/oj). 



small energy E. In this limits SW orbits are elliptical 
and the integrals in Eqs. (I12I) may be easily performed, 
leading to 

Sl[E) - fl (1 - E/E ) ; (17) 
U(E) = uE; V{E)=vVE; W{E) = w/Ve, 

where 

ft = jH z VT+h, E =E b (4 + 4h)/(2^ 

v = ^\f~E h \f\ 
~h) /h. 



h), 



u = fi (2 + h)/Vl + h, 

7/^ (l + h-VT 



and h = H x / H z . We have checked numerically that these 
expressions remain almost quantitatively correct in the 
entire relevant range of parameters. With these forces 
one finds J{E) — E a , where a = (2w — v)/2v is eccen- 
tricity index which ranges from a — in the absence of 
the easy plane anisotropy h — > to a = —1/2 in the 
strong easy plane case h >• 1. 

With these approximations the top of the separatrix 
line, given by the conditions % — and <j) — 7r, is found 
from the cubic equation for E 1 ^ 2 , cf. Eq. (ITU) . 
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E 3/2 = wl s . (18) 



It may have either three or one solution with the bifurca- 
tion point determining the LU2(I S ) relation. The critical 
current (fTl?)) is then found by the substitution of E = E^ 
in the left hand side of Eq. (|18|) . This way one finds an 
approximate location of the two lines depicted in Fig. |H1 
Their intersection results in the theoretical lower bound- 
ary for swing-switching critical current, which happens 
to be practically ^.-independent for h > 1 



X s >H z /24. 



(19) 
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FIG. 10: (Color online) Energy E vs time (in units of 
(■yMs)' 1 ) for Q = 100 (blue, dot-dashed), Q = 2 (black, 
solid), and Q = 1 (red, dashed). Same parameters as in Fig[2] 
except of a, which is changed to adjust Q. 

One may notice that in our simulations the critical 
switching current is about factor of 3 larger than this 
idealized estimate. The reason is that the estimate based 
on the conserved quantity (fTil) neglects the dissipation. 
Calculating the work of the friction force as described 
below Eq. ([HI), one finds E d cx E b a(2 + h)/y/l + h with 
the proportionality coefficient which logarithmically di- 
verges at uj — > u>2 (indeed, at the bifurcation point the 
dissipation always prevents the switching). In practice 
we found from the simulations that the condition 

e = ^^±5>i ( 2 „) 

ba 2 + h 

provides a good estimate for the underdamped range 
(Ed/ Eb < 1). If Q ^> 1 one expects to approach the 
theoretical boundary (fl9|) . In the opposite case Q < 1 
the AC perpendicular method losses its advantage over 
the DC approach, because the effective oscillator is over- 
damped. 

One may notice that the simulations discussed above 
happen to have a mediocre quality factor Q w 2, which 
may explain why they result in the critical current which 
is bigger than the theoretical bound (TT9|) . To test these 
predictions we performed additional simulations. Fig- 
ure [TO] shows E(t) for Q = 1, 2, and 100. Notice that 
for Q — 1 E(t) is overdamped and for Q = 100 E(t) 
is highly underdamped. Notice also that in the under- 
damped Q — 100 case the maximum energy E max is 
significantly higher than in the other two. Expanding 
this analysis further Fig. [TT] shows the maximum tran- 
sient energy E max as a function of h for several values 
of Q. To keep Q constant we adjusted dissipation a ac- 
cording to Eq. PU|) . while H z and I s arc kept fixed and 
u) = Qo(h). Notice that the maximal energy remains al- 
most completely constant for each value of Q and grows 
with increasing Q, so that for Q > 1 the maximum en- 
ergy is nearly 50% larger than for Q = 2. This shows 
that, if the spin-current is scaled with H z and AC fre- 



10 20 30 40 50 60 
h 

FIG. 11: Maximum transient energy E ma x/Et as a function 
of h = H x /H z for several values of Q. Here H z = 0.033M 3 , 
Is = 0.003M fl , T = 0, and u = Q - 



quency is scaled with fig , the system behavior is governed 
by the single scaling parameter Q, defined in Eq. (|20p . 
This observation reduces the multidimensional space of 
system's parameters down to the single essential param- 
eter: the quality factor Q. 



IV. CONCLUSION 

We have shown that by using an AC spin current with 
driving frequency close to the precession frequency of the 
free layer one may achieve the reliable switching of spin- 
torque valves. Kee ping advantages of the DC perpendic- 
ular approach [l3|, Il4 l2ll . |22| , the switching times are 
well below Ins. In addition, the AC method offers a sig- 
nificant reduction of the critical switching current and 
associated energy dissipation. We explained the effect 
by the underdamped nature of the effective non-linear 
oscillator, which allows to cross the energy barrier by 
transient oscillations. The latter may significantly over- 
shoot the equilibrium energy, reducing the spin-torque 
required for the switch to occur. 

We have developed a theoretical description of the ef- 
fect which allowed us to identify the optimal AC fre- 
quency with the upper bifurcation frequency of the driven 
non-linear oscillator. It also provides with the estimate 
of the current amplitude needed to perform the switch. 
Possibly most essentially, it provides the scaling argu- 
ments showing that I s cx H z ; uj m Q — -fH z ^/l + h, 
while the relative advantage of the AC method is gov- 
erned by the single scaling parameter: the quality factor 
Q, Eq. O. 
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